超临界环境下煤油替代物液滴燃烧特性数值研究

为研究液氧/煤油火箭发动机燃烧室内经喷注形成的煤油液滴的燃烧过程，基于实际气体状态方程、高压热物性修正、高压气液平衡和详细化学反应动力学，建立一维的全瞬态液滴燃烧模型，对超临界环境下两组分煤油替代物液滴的燃烧特性及液滴初始直径的影响进行仿真研究。结果表明，在超临界环境下，相比于煤油液滴纯蒸发过程，煤油液滴燃烧过程的迁移时刻大大提前；煤油液滴着火之后很快进入超临界燃烧阶段，此时液滴燃烧过程可以看成中心附近的燃料高浓度区与外侧氧气高浓度区之间的扩散燃烧过程；煤油液滴的火焰半径先增大，达到最大值之后开始减小，并减小为零，火焰温度在着火之后快速上升至最大值，并基本保持不变，在火焰半径减小为零之后开始降低；随着液滴初始直径的增大，火焰特性以及液滴中心参数变化曲线趋势不变、整体延迟，着火时间、迁移时间和液滴寿命增大。     

类X-37B升力体再入飞行通信特性中断研究

针对再入非平衡等离子流场的数值仿真方法及等离子体鞘套对不同波段电磁波传输特性的预测开展研究。首先，利用某工程钝头体飞行器再入实测通信飞行测量数据，初步校验数值仿真等离子体特性和通信特性方法的正确性；其次，在典型弹道下针对升力体再入飞行器开展等离子体鞘套数值仿真，分别给出S波段、Ka波段与K波段电磁波衰减特性，研究指出结合当前国内中继星体制情况采用高频段电磁波Ka波段或K波段中继通信方式是解决通信黑障的技术途径之一；最后，研究了电磁波数衰减受传播方向角和攻角的影响规律。研究结果可为发展解决我国再入飞行器通信黑障的设计方案提供参考。     

非预混条件下的旋转爆轰燃烧室双波头演化过程数值模拟

针对旋转爆轰燃烧室双波头演化过程中流场结构变化的问题,对非预混条件下的旋转爆轰燃烧室从起爆到形成稳定的双波头过程进行了数值模拟研究。研究结果表明,从起爆到形成稳定爆轰过程,燃烧室主要经历了起爆、爆轰波对撞和稳定爆轰三个阶段；在爆轰波对撞阶段,首次对撞是两个爆轰波间的对撞,由于对撞点处缺少新鲜混合气,从而在对撞结束后衰减为两个压力波。第二次对撞是两个压力波间的对撞,因为在第二次对撞点附近存在新鲜混合气来支撑爆轰波的持续传播,故对撞结束后产生了一个爆轰波和一个较弱的压力波；第二次对撞发生后,燃烧室内的压力波反射叠加并形成局部高压区,此高压区压缩气体使气体温度升高,高温气体引燃混合气后,最终发展成为第二个爆轰波；稳定阶段,两个爆轰波均能稳定自持传播,爆轰波峰面压力可达1.45MPa,波后温度为2500K,爆轰波速度稳定在1738m/s,产生的推力与比冲分别为79.76N和2312.15s；斜激波的存在使燃烧室出口平面流场产生了较大波动。      

A comprehensive study of the 2016 Mw 6.0 Italy earthquake based on high-rate (10 Hz) GPS data

The 2016 M_w 6.0 Italy earthquake is successfully recorded by the near-field 10?Hz GPS and 200?Hz Strong Motion (SM) stations, providing valuable data for this study. A comprehensive study of this earthquake is carried out based on GPS data, which contains coseismic deformations analysis, noise analysis, seismic wave picking, and magnitude determination. The noise of most GPS-derived displacement waveforms can be described as a combination of white noise, flicker noise, and random walk noise after the earthquake occurrence, and the spectral indices vary significantly for most stations, implying that the seismic signals have affected the noise characteristic of GPS-derived displacement waveforms. S-transform is employed to assess the GPS capability to detect the seismic arrival time. The SM station AMT and the GPS station AMAT are in good agreement in seismic wave picking, and the difference is only 1.2?s in the north component, suggesting that the outcome of seismic wave picking using GPS data is reliable. Then, a classic empirical formula is employed to determine the moment magnitude. A robust moment magnitude (Mw 5.90) can be estimated by the nine GPS stations with about 23.9?s. If four GPS stations near the epicenter is chosen to determine the magnitude, it only take 13.0?s to retrieve a reliable preliminary (Mw 5.82) magnitude, which is 5.4?s ahead of nine stations. In addition, Cross Wavelet Transform (XWT) is adopted to measuring the correlation and phase relationship between GPS and SM records. The result of XWT analysis indicates 10?Hz GPS is capable of capturing reliable and accurate coseismic dynamic deformations, as evidenced by the XWT-based semblance being close to 1 between GPS and SM records. The above results confirm the capability of 10?Hz GPS to capture coseismic dynamic deformations, detect seismic arrival time, and determine earthquake magnitude. Moreover, rapid magnitude determination based on 10?Hz GPS data can be regarded as an important supplement to Earthquake Early Warning (EEW).     

预旋喷嘴径向角度对预旋特性影响的数值研究

为了降低低位进气预旋流路的气动损失，针对带有不同径向角度（0°~30°）预旋喷嘴的预旋系统进行了数值仿真，并对流动特性、温降特性和比熵增特性进行了分析。结果表明:随着预旋喷嘴径向角度的增大，预旋系统无量纲温降先增大后减小，流动阻力减小，预旋系统的流量随之增大。旋转雷诺数为2.3×107时，预旋喷嘴带径向角度的预旋系统无量纲温降比传统喷嘴最大可提高18.3%，存在某一角度使预旋温降特性达到最好。预旋系统内的耗散主要发生在预旋腔和共转盘腔内，径向角度为10°时其比熵增变化量分别占整个预旋系统总体比熵增的42.4%和30.2%；合理设计预旋喷嘴的径向角度，能改善预旋腔内气流的流动效果，并且可以减少整个预旋系统的不可逆损失。      

基于变刚度变阻尼的柔性喷管动态模型

针对柔性接头动态迟滞曲线受控制系统控制位置精度和动态响应速度影响较大的问题,基于电液伺服机构和柔性接头变刚度变阻尼模型,构建了柔性喷管的电液伺服机构-变刚度变阻尼模型,将其和电液伺服机构-定刚度定阻尼模型进行了对比。分析了电液伺服机构主要参数、柔性接头工作参数等对电液伺服机构-柔性接头系统动态特性的影响。分析结果表明:电液伺服机构-变刚度变阻尼模型所构造的迟滞曲线可更准确地与实验结果相吻合,并符合迟滞曲线随频率变化的规律,反馈系数、放大器静态放大系数、电液伺服机构增益、滑阀流量增益等参数对系统动态特性的影响更为明显。该模型为固体火箭发动机电液伺服机构-柔性接头系统动态特性的调整提供理论依据。     

微阴极电弧推力器羽流探针诊断研究

微阴极电弧推力器（micro-cathode arc thruster,μCAT）具备功率低和结构简单的特点,能够满足微纳卫星的任务需求,具有良好的发展前景。μCAT羽流的诊断可以揭示推力器的加速机理,对提高其性能具有重要意义。利用朗缪尔三探针对μCAT羽流进行诊断,得到了μCAT羽流不同位置的电子温度、电子密度和离子速度等羽流特性,研究了外加磁场、充电时间和阴极材料对羽流特性的影响。研究结果表明,μCAT放电初期产生的等离子体电子温度较高,密度较大;随着等离子体向下游运动,电子温度和电子密度降低,离子速度增大;外加磁场的磁感应强度越强,电子温度和离子速度越高,电子密度有所降低;磁场位置适当向推力器下游平移,能够有效提高推力器中轴线的电子密度;μCAT充电时间越长,电子温度、电子密度和离子速度越大;相比于CuW和AgW阴极,Ti阴极羽流的电子温度更高,电子密度更低。     

平板式预膜喷嘴初次雾化特性试验

为深入了解平板式预膜喷嘴初次雾化特性,试验研究不同进口条件对平板式预膜喷嘴初次雾化特性的影响规律。试验采用了背光照明和片状激光照明相结合的高频拍摄手段分别获得液膜俯视和侧视破碎形态,同时运用本征正交分解(POD)法和液膜边缘定位等分析方法进行光学图像结果后处理,获得表征平板式预膜喷嘴初次雾化特性三个物理量:液膜波动频率、破碎距离和横向不稳定波长。试验结果表明:①通过分析侧视和俯视破碎形态,平板式预膜喷嘴液膜破碎形态可分为三类:末端破碎、波浪脱落和表面剥离,进口韦伯数对预膜喷嘴破碎形态的影响占主导地位;②把POD法和液膜边缘定位方法等相结合方法应用到高频非接触光学喷嘴雾化图像的后处理分析中,是一种非常有效的数据后处理方法;③液膜初次雾化特性主要受到进口韦伯数和气液动量比的影响,破碎距离和横向不稳定波长都随进口韦伯数的增加而降低,液膜表面不稳定波动频率随进口韦伯数的增加而增加,所获得的经验关系式与试验数据吻合较好。所获得的进口参数对液膜破碎形态和雾化特性的影响规律为喷嘴后续优化设计提供了依据。     

推力室状态对其声学振型及其阻尼特性的影响

为确定工程中冷态条件下获得的推力室声学特性能否表征真实条件下的声学特性,研究了冷态无流动、热态气相流动和湍流两相燃烧三种状态下推力室声学振型及其阻尼特性。在推力室稳态流场中的有限区域施加数值定容弹,激发其具有多模态声学振型的大幅值压力振荡,采用衰减时间和半带宽来定量评价所激发的不同声学振型压力振荡衰减快慢,进而获得其阻尼特性。在相同过载比的数值定容弹激励下,在冷态条件下能激发包含更多声学振型压力振荡,且该振荡衰减时间更长,相同振型压力振荡衰减比热态条件下慢。在冷态条件下,一阶切向振型振幅最大,为最容易被激发声学振型;一阶纵向振型半带宽最小,为最难衰减的振型。在热态条件下,一阶纵向振型为最容易激发声学振型,也为最难衰减声学振型。从所激发的主要振型及其相对衰减的快慢来看,冷态条件下获得的声学特性能够表征真实条件下的推力室的声学特征。     

纳秒脉冲表面介质阻挡放电激励改善飞翼力矩特性的实验研究

针对飞翼布局力矩控制问题,采用纳秒脉冲表面介质阻挡放电（NS-DBD）激励,在来流风速30 m/s时,开展飞翼等离子体流动控制风洞试验,研究了不同激励参数和位置对飞翼升阻特性和力矩特性的影响。结果表明,NS-DBD激励能够有效改善飞翼大迎角气动特性。激励频率对飞翼升阻特性影响较大,激励频率为0.2 kHz时,增升效果最好,最大升力系数提高14.5%,失速迎角推迟5°。随着激励频率的增加,增升效果逐渐变差,减阻效果变好。单侧施加激励时,能够实现大迎角下飞翼模型的力矩控制,随着激励频率的增加,滚转力矩的控制效果减小,激励频率为0.2kHz时,平均滚转力矩系数变化为ΔMX=0.005691;偏航力矩的控制效果增大,激励频率为1kHz时,平均偏航力矩系数变化为ΔMY=-0.001571;俯仰力矩的控制效果减小,激励频率为0.2kHz时,平均俯仰力矩系数变化为ΔMZ=-0.002576。在中翼段和内翼段施加激励,破坏了飞翼的俯仰力矩特性,在外翼段和机翼右侧施加激励,能够显著改善飞翼的俯仰力矩特性。流场测量结果表明：等离子体激励对飞翼气动力矩的控制,主要是通过控制流动分离和控制横向流动来实现的。NS-DBD激励为改善飞翼布局稳定性和操纵性提供一种潜在的技术手段。